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Abstract — For a right-invariant and controllable driftless sys- 
tem on SU(n), we consider a time-periodic reference trajectory 
along which the linearized control system generates 5u{n): such 
trajectories always exist and constitute the basic ingredient 
of Coron's Return Method. The open-loop controls that we 
propose, which rely on a left-invariant tracking error dynamics 
and on a fidelity-like Lyapunov function, are determined from 
a finite number of left-translations of the tracking error and 
they assure global asymptotic convergence towards the periodic 
reference trajectory. The role of these translations is to avoid 
being trapped in the critical region of this Lyapunov-like 
function. The convergence proof relies on a periodic version 
of LaSalle's invariance principle and the control values are 
determined by numerical integration of the dynamics of the 
system. Simulations illustrate the obtained controls for n = 4 
and the generation of the C-NOT quantum gate. 

I. INTRODUCTION 
Consider the right-invariant driftless system 

m 

X = Y,UkHkX, X{0)^I, (1) 
fe=i 

where X £ M" is the state, A/" is the Banach space 
of square n x n matrices with complex entries endowed 
with the Euclidean norm, H = {Hi,...,Hm} C su(n), 
M^. e M are the controls, and / is the identity matrix of 
M". The periodic motion planning problem for this system 
is formulated as follows. Given a goal state Xoo G SU(n) 
and T > 0, find a smooth periodic reference trajectory 
Xr. M+ SU(n) of period T, with ^^(0) = X^o, and 
determine continuous open-loop controls ut'- M+ M, for 
1 < /c < m, in a manner that the tracking error between the 
trajectory X: M+ SU(n) of ([T]) and Xr converges to zero 
as < ^ cx), that is, \mit^oo[X{t) — X^it)] = 0. 

We remark that there is no loss of generality in as- 
suming that X{Q) = / in ([T]). Indeed, since system 
([T]l is right-invariant, if (ui(t), . . . , for 

t e M+, is a solution of dUi with X{0) = /, then 
{X{t)Xa, {ui{t), . . .,Ujnit))), for t e M+, is a solution of 
^ with initial condition X{0) = Xq G SU(n). Therefore, 
if the periodic motion planning problem has been solved for 

The first author was fully supported by CAPES. The third author was 
partially supported by CNPq. The second and third authors were partially 
supported by CAPES/COFECUB and "Agence Nationale de la Recherche" 
(ANR), Projet Blanc CQUID number 06-3-13957. 

H. B. Silveira and P. S. Pereira da Silva are with Laboratory of Automa- 
tion and Control, Department of Telecommunications and Control Engineer- 
ing, University of Sao Paulo, Brazil hectorbessaSyahoo . com . br 
and paulo@lac.usp.br 

P. Rouchon is with Mines ParisTech. Centre Automatique et Systemes, 
Mathematiques et Systemes, France 
pierre . rouchon@mines-paristech . f r 



system ([T]i with X(Q) = I, it is straightforward to show that 
it will also be solved for ^ with X{0) = Xq e SU(n). 

The main result of this paper is the determination of 
a solution for the periodic motion planning problem. This 
is established by Theorem |2] in Section 2, whose only 
assumption is that system ([U regular, in the sense of 
Definition [T] in Section 2. The results of Coron's Return 
Method show that such condition is always met in case 
the system is controllable on SU(n) (see Remark |2] in 
Section 2). Loosely speaking, by finding an appropriate 
reference trajectory Xr, using the time-dependent change of 
coordinates Z = Z{X,t) = X'^Xr{t), which corresponds 
to the tracking error on the group SU(n), and defining an 
adequate "feedback", we determine an algorithm that obtains, 
in a finite number of steps, continuous open-loops controls 
Uk, for every 1 < A; < m, which assure that the tracking 
error X — Xr converges to zero as < ^ oo. This algorithm 
relies on Lyapunov-like convergence results inspired in the 
periodic version of LaSalle's invariance principle presented 
in [11], and in the ad-condition stabilization method of [6]. In 
a certain sense, we have used the real part of the trace of the 
left-invariant tracking error Z as a Lyapunov-like function, 
that is, V{Z) = ^{tr{Z)). In the case of quantum systems, 
V can then be seen as a fidelity-like Lyapunov function. 

The problem of steering a quantum system from a given 
initial state to an arbitrary final state, which can be regarded 
as a particular case of the periodic motion planning problem 
here formulated, has recently been treated in [8] using a 
flatness-based approach and in the book [4] (see also the 
references therein), where many quantum control techniques 
used in the literature are grouped together and explained in 
detail, such as Lyapunov-based methods, optimal control and 
decompositions of SU(n). Our Lyapunov-like approach has 
no restrictions on the goal state Xo^ & SU(n) and on n, as 
long as system (HJ is regular 

The layout of the paper is as follows. Section 2 is 
entirely dedicated to the proof of Theorem |2] mentioned 
above. Simulations illustrate in Section 3 the generation of 
the Controlled-NOT (C-NOT) gate for a quantum system 
with n = 4. Appendix presents the proof of the important 
convergence result of Theorem [T] in Section 2. 

II. Main Result 
Based on ([T]i, we define the reference system 

m 

Xr^Y. <HkXr, Xr{0) ^ X^ G SU(n), (2) 

k=l 



where Xr G Af" and the smooth time functions u^^,: M ^ M 
are still to be specified. 

Definition 1: System ([T]) is said to be regular when, given 
T > 0, there exist smooth periodic functions u^: R ^ M 
of period T, for all 1 < fc < m, such that the solution 
Xj: R+ ^ SU(n) of ©, with Xj{0) = / and = u^, 
is also periodic of period T and satisfies 

span{B^(0),l < fc < m, j e N} ==su(n), (3) 

where N is the set of natural numbers (including zero), 

Bi+\t) = -A{t)Bl{t)+Bi{t), 1 < fc < m, J G N, t e R. 

Remark 1: Note that A: R ^ 5u(n), B^, B^,: R M" 
are smooth and also have period T, for every 1 < k < m, 
j G N. Hence, they are bounded mappings. 

Remark 2: Note that the hnearized control system of (|2|i 
(or of ([T]i) along the trajectory {X^ , {uj , . . . , u^)) is given 
by = A{t)X^^ + jZlLiWkBlit), Wk £ R. Based on 
Coron's Return Method (see [2], [3]), it can be shown that 
([TJ is regular in case Lie(iJ) = su(n). We recall that ([TJ is 
controllable on SU(ri) if and only if Lie(iJ) — su{n) [1]. 

For simplicity, we shall assume throughout this paper that 
system ([TJ is regular, that T > has been fixed and that 
the functions u^, in dU were specified accordingly, that is, 
= uJI, for 1 < fc < TO. Moreover, we also assume that 
the goal state X^o e SU(n) is fixed. Define Xr'. K ^ SU(n) 
as Xr = X^Xoo- Note that Xr is the solution of (|2|i 
with Xr{0) = Xoo and that Xr also has period T. It 
will be shown afterwards that Xr can indeed be used as a 
reference trajectory. We also adopt the following notations. 
The imaginary unit of C is denoted by i and if z G C, then 
^{z) is its real part and ^(z) its imaginary part. 

It is straightforward to verify from ([T]i and (|2]i that the 
time-dependent change of coordinates 

Z = Z{t,X) = X^Xrit), for all {t,X) G R x M", 

along with the time-varying control shift 

Vk = u^{t)—Uk = u^{t)—Uk, for all i G R, 1 < fc < m, 

determine the left-invariant "closed-loop system" 

Z = ZXlit) VkHkXrit), Z(0) - G SU(n), (4) 

k=l 

for all {t, Z) G R X M". If we can find continuous functions 
Vk- R+ R, for each 1 < fc < m, such that 

lim Z{t) = lim X''{t)Xr{t) = /, (5) 

t — >oo t — ^oo 

where Z: R+ SU(n) is the solution of system (|4|l and 
X: R+ SU(n) is the solution of system ([TJ with the 
continuous open-loop controls 

Uk{t) = ul{t) - Vk{t), for all t G R+, 1 < fc < to, 

it is then clear that 

lim [X{t) - Xr{t)] = 0, (6) 

t — >oo 

thus solving the periodic motion planning problem. 



Let V: Ar ^ R be defined by 

V{X) ^ 5R(ir(X)), for all X G M", (7) 
and consider the auxiliar system 

W - WXl{t) J2 fkdkit, W)HkXr{t), (8) 
fc=i 

where {t, W) G R x A/", ^ is a fixed real number, 
1 < fc < m, and 

ak{t, W) = fkXiWXl{t)HkXr{t)). (9) 

Notice that the "closed-loop" system ((4j with "feedbacks" 
Vk = fkCikit, Z) is nothing but the auxiliar system ((8j-((9j. 
Note also that in ((TJ is linear- and that, for X G SU(n), 
we have —n < V{X) < n and V{X) = n if and 
only if X = /. Furthermore, by construction, V{t, W) = 
EfeLi W)^ > 0, for all {t, VF) G R x M". 

In what follows, we shall show how the next theorem, 
which is a Lyapunov-like convergence result for the auxiliar 
system with Lyapunov-like function — ^{tr{W)), 

and whose proof is deferred to Appendix, determines con- 
tinuous functions Vk '■ R+ R, for 1 < fc < to, such that ((5j 
is satisfied for the "closed-loop" system ((4j. We remark that 
the properties of V stated above are essential in the proof. 
Our approach to solve the periodic motion planning problem 
is then summarized in Theorem [2] 

Theorem 1: Consider the set 

n 

G = {a;GR:a; = ^K(A,), for some Ai G C such 

i=l 
n 

that |A,| = 1, J] A, = 1, 3(Ai) = • • • = 3(A„)}. 

1=1 

Then, G is a finite set, n E G and n — max(G). Further- 
more, letting 6 be the maximal element of the set G \ {n}, 
we have that, for all q = (to, WtJ G R x SU(n), 

where Wq. R -» SU(n) is the solution of ([8j-((9j with initial 
condition W^g(io) — Wt^- 

Suppose that t/(Xoo) > 5- In Theorem [T] we choose 
q = (0,Xoo) G R x SU(n). Therefore, Ymit^oo Wq{t) = I. 
Hence, the smooth "feedbacks" Vk'. R+ R defined as 

Vk{t) = fkak{t,Z{t)) - f^Y{xHt)HkXr{t)), 

for i G R+, 1 < fc < m, assure that Z{t) = Wq(t), for 
< G R+. Indeed, compare (gj with ([8j-((9j. Thus, ((5j holds. 

Now, assume that V^(Xoo) < S. For this case, based 
on continuity arguments, we determine an adequate (con- 
tinuous) path from Xoo to I which, in a certain sense, 
reduces the problem to the situation where y(Xoo) > S. 
In order to achieve this, the main idea is to find a path 
Z: [0,1] ^ SU(n), with Z{0) = Xoo and Z(l) = /, 
and obtain ^ Oq < 9i < ■ ■ ■ < O^-i < 9^ ^ 1, 
such that V(Z{di+i)'<'Ziei)) > ^, for all < £ < X - 1. 
It thus follows from Theorem 1 that, for 1 < i < N, 



limt^oo Z{9i)Weit) = ^(6*^), where Wf. M ^ SU(n) 
is the solution of (O-® with initial condition We{Te) — 
Z{ei)^Z{0i_i) G SU(n), where = Ti < • • • < Tn+i are 
such that We{Ti4-i) w /. Loosely speaking, we then "glue" 
together the left-translations ^(^1)1^1, . . . , Z{9n)Wn in an 
appropriate manner in order to define a continuous solution 
{Z(t),{vi{t),...,Vm{t))), for t e M+, of system dUl that 
satisfies (|5]l. We remark that, for every 1 < < A^, it is as if 
we were in the case V{Xoo) > S. In the sequel, we formalise 
these arguments in detail and determine an algorithm which 
obtains, in N steps, continuous functions v^: R_|_ R, for 
1 < k < m, such that ^ holds. 

It is a standard result that any Xoo G SU(ri) can be written 
as Xoo — Mtdiag(expiAi, . . . ,expiA„)M, where M is a 
unitary matrix, Ai, . . . , A„ G M and — ^- Consider 

the path Z: [0, 1] SU(n) from Xoo to / defined by 

Z{e) = A/tdiag(expiAi(l ~9),.. . ,expjA„(l - e))M, 

for all e e [0,1]. Let a,b e [0,1]. Hence, Z{byZ{a) = 
M^diag(exp «Ai (5 — a), ... , exp iXn{b — a))M and there- 
fore V(Zib)1Z{a)) = E"=iCOs(Aj(fe - a)). Since the 
function 7: [0, 1] ^ R defined by 7(0) = cos(Aj6'), 
for all 9 £ [0, 1], is continuous with 7(0) — n, there exists 
> such that -f{9) > 6 in case [^l < i^, for all 9 E [0, 1] 
(indeed, choose e = n - (5 > 0). Hence, V(Z{by^{a)) > S 
whenever |6 — a| < ly, for all a, 6 e [0, 1], and there exists a 
non-zero r/ G N such that, for all N > rj, 

n 

V{z\+;Zf) = cos(A, A) > S, (10) 
i=i 

for all < £ < TV - 1, where Ze = Z{9e), 9e ^^A, 
for every < ^ < iV, with A = l/N. Note that Zq = 
Z{0) = Xoo and Zn = Z{1) /. Let iV > 77 and consider 
the continuous function f3: Af" x A/" R defined by 
l3iX,Y) = y(ytx), for all iX,Y) G Af" x A/". Since 
SU(n) xSU(7i) is compact, /3|(SU(n) xSU(n)) is uniformly 
continuous. Therefore, by ( fTOb . there exists /i > such that, 
for all X G SU(n) and < ^ < iV - 1, we have 

\\x -Z4< ^i^v(zl+^x)> 6 (11) 

(indeed, choose e = cos(AjA) — (5 > and consider 

the sup norm on Af" x Af"). The aforementioned algorithm 
is described below. Recall that Zq = Xoo and Zn ~ I- 

Algorithm 1: Let Xoo G SU(n). Choose any non-zero 
iV G N in a manner that ([Toll holds. Define Ti = and 
Wo(ri) = /. For every I < £ < N — 1, choose a real number 
Te+i > Ti such that V (zl^,^ eWi{Te+i)) > S, where 
Wf. R — + SU(n) is the solution of the auxiliar system (O-® 
with initial condition W^(rf) ='zl'Ze^iWe-i{Te) G SU(n). 
Define 

Z{t) = ZiWi{t) G SU(n), for t G [rf,Tf+i), 

i^fc(t) = fkak{t, Wi{t)) G R, 1 < fe < m, for t G [T,, T.+i). 



If Tat > Tat-i has been chosen as above, define 

Z{t) = WN{t) eSU{n), for t>TN, 

Vk{t) = fkakit, WNit))eR, 1 < fc < m, for i > T^, 

where Wn- R SU(n) is the solution of dSj-® with initial 
condition WNiTN) = Zn-iWn-i{Tn) e SU(n). ■ 
Some remarks are in order First of all, from the rea- 
soning preceding Algorithm 1, we know that there al- 
ways exists some non-zero iV G N such that (fTOl i is 
true. Furthermore, Theorem [T] and property ( fTTT i assure 
that T^+i > Tf can always be chosen as required in the 
algorithm, for every 1 < £ < A^ — 1. It is also clear 
that {Z(t), ivi{t), . . . , Vm{t))), for t G R+, determined by 
the algorithm is a continuous solution of the "closed-loop" 
system (|4]l. Indeed, compare (|4]l with (l8]l-(|9]l. Finally, since 

V(Z^j^ZN-iWN-iiTN)) = V(ZN-iWN-iiTN)) > S, 
Theorem [T] implies that \imt^ooW]y{t) = I. However, 
Z{t) = WN{t), for t > Tjv. Therefore, the continuous 
functions Vk '■ R+ — > R determined by Algorithm 1 are such 
that (|5]l is satisfied. We have thus shown our main result: 

Theorem 2: Assume that system ([T]) is regular, in the sense 
of Definition [T] Given Xoo e SU(n), T > and "feedback 
gains" fl > 0, consider Xj: R SU{n) and u^: R ^ R as 
in Definition[Tl for 1 < k < m. Define X^ ~ Xj Xoo- Then, 
there exist continuous open-loop controls uj,: R+ — * R, 
for 1 < k < m, such that ^ is satisfied. In other 
words, the periodic motion planning problem always has a 
solution when ^ is regular More precisely, if V{Xoo) > S, 
where S is as in Theorem [T] then the smooth open-loop 
controls Uk{t) = ul{t) - f^Y(X'^ {t)HkXr{t)), obtained 
by numerical integration, for t G R+, 1 < k < m, assure 
that (|6]l holds. Otherwise, in case V{Xoo) < S, then by 
following Algorithm 1 we determine continuous functions 
Vk- R+ R, for 1 < fc < m, such that the corresponding 
continuous open-loop controls Uk{t) = uj^{t) — Vk{t), for 
t G R+, assure that ^ is satisfied. 

III. QUANTUM MECHANICAL EXAMPLE 

After some approximations, an appropriate change of 
coordinates, scalings and simplifications, a controlled quan- 
tum system consisting of two coupled spin-^ particles with 
Heisenberg interaction and driven by an external electromag- 
netic field, can be modeled as [4] 

Y ^{D + D,u, + DyUy + D,u,)Y, Y{0) ^ I, (12) 

where Y G M'^ (n = 4), the controls u^., Uj,, G R are the 
X, y, z components of the electromagnetic field, respectively, 
D = diag(3i, -I, ~t, -i), = - 3i7j|, Dy = + 
3Hi^, D, = iJ/l - 3i?ii G su(4), and H^^ = (h^f), 
HJj — {hl:^-') G su(4) are the matrices with entries 

hf/' = 1, - -1, = 0, for k,i^i,j, 

hlf = h^f = I: hlf = 0, for kj^i, j, 

respectively, for all 1 < i < j < n. 

Now, in order to remove the drift term DY in ( fT2] i. we 
define, as usual, the time-dependent change of coordinates 



X = $(t, Y) = e--°*r, for alWi, Y) e 



[I a, : 
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X M^. In these 



coordinates, ( fT2] i is described 

X = (C^Uj; + CyUy + C X , 

where 



X(0) 
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for all t G M. We choose the real controls u^, Uy, as 





= (ui - 
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respectively, for all t e R, where ui, . . . ,uq G M are the 
new controls. By applying the rotating wave approximation 
(RWA) (see e.g. [9], [4], [5]) to system (fT3])-(fT4b. which con- 
sists in considering only the terms that are time-independent 
and in disregarding all the oscillating ones, we obtain the 
following time-independent driftless system 

(15) 

with initial condition X{Q) = /. It is straightforward to ver- 
ify that Ut{{H?;,.Hi^,H^^,Hi^,H^^,Hi^}) = su(4), i.e. 
the system is controllable on SU(4). Hence, Coron's Return 
Method implies that the system is regular (see Remark [21) 
and therefore Theorem |2] can be applied. We choose T = 1 
and as goal state the C-NOT (ControUed-Not) gate 



/ 1 



V 





-1 

/ 



e SU(4), 



which is one of the universal gates and has great importance 
in quantum information theory [7], [5]. It is easy to see from 
the proof of Theorem [T] that G = {-4,0,4} with 5 = Q. 
Since V{Xoo) — 2>0, Theorem |2] implies that the smooth 
open-loop controls uu{t) = ul{t) - f^Y{X^t)HkXrit)) = 
u\{t) — for t G R+, 1 < fc < 6, obtained by numerical 
integration, assure that \imt^oo[X{t) — Xr{t)] = 0, for any 
"feedback gains" > 0. Here, uj. = itf with T = 1, 
and Hi 
H5 



ttR 
-"12' 



ttR 
^li' 

Ha 
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'■k 
H? 



3 - -^13, i?4 

However, the periodic functions 



-f^l3' 



In quantum mechanics, this description is usually called the interaction 
picture or interaction representation. 



u\, . . . ,u\ are not known explicitly. Coron's Return Method 
only establishes their existence. Fortunately, for system ( fTSl ). 
symbolic computation software packages have shown that if 
we define them as u\{t) — X^^ii afc£ sin(27r£i), for i G M, 
1 < A: < 6, with nj > \ and where aki G M are randomly 
chosen from the uniform distribution on the interval [—a, a] 
with "sufficiently large" a > 0, then it is "very likely" that 
dim(span{S^(0), 1 < A: < 6, < j < 6}) = 15, that is, © 
holds (recall that dim(su(4)) = 15). And, when (O is true, 
it follows that \iiTLt-^oo[X{t) - Xr{t)], where Xr = X^Xoo- 
We remark that since u\ is an odd periodic function with 
period T = 1, the solution X^.: M SU(n) in Definition[T] 
is also periodic with period T = 1. Note that a and n/ 
determine the "excitation level" of u].. For fk = 1, computer 
simulations have suggested that as a and Uf get larger, the 
faster the convergence of the tracking error X ~ Xr to zero 
(assuming that dim(span{i3*^(0)}) = 15, of course). 

The obtained simulation results are now presented for 
/fc = l,a = 7i/ = 5 and aki having as values the 
corresponding entries of the matrix A = {ake) below 
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\ -0.18 



-1.39 
-1.54 
-0.31 
0.77 
0.60 
-4.44 



4.66 
0.92 
1.75 
4.09 
-0.27 
-1.38 



4.31 
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-3.75 

2.59 ) 

verified 



that 



With these choices, we have 
dim(span{_B|'(0)}) = 15. Figure 1 exhibits the convergence 
of \X — Xj\ to zero (Euclidean norm on M**). We see that 
the norm of the tracking error is non-increasing. In Figure 2, 
the controls ui, U2 (top) and the "feedbacks" v\,V2 (bottom) 
on the time interval [0, 10] are shown. Notice that is 
relatively small in comparison with the control Uk, for 
fc = 1,2. Therefore, the control Ufc is relatively close to 
as defined above, for fc = 1,2. In order not to overwhelm 
the presentation, we have chosen not to exhibit Uk,Vk, for 
3 < fc < 6. They have, however, a similar behavior and a 
similar order of magnitude as for fc ~ 1,2. 
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Fig. 1. Convergence of the norm of the tracking error to zero. 

IV. CONCLUDING REMARKS 

In the solution here presented for the periodic motion 
planning problem, the only needed assumption is that system 




'oi 2345678910 

t 



Fig. 2. Controls ui, J12 and "feedbacks" vi,V2 on the interval [0, 10]. 

([T]i is regular, which requires that the periodic functions 
satisfying ^ are explicitly known. Nevertheless, this will 
hardly be the case in general. For this reason, currently under 
investigation is the explicit determination of u^^ in (|2]l in a 
manner that Theorem [T] still holds under assumptions other 
than the regularity of system ([T]). 

V. ACKNOWLEDGMENTS 

The authors would like to thank Jean-Michel Coron and 
Mazyar Mirrahimi for valuable discussions and suggestions. 

Appendix 

In order to prove Theorem [T] we need first a few inter- 
mediate definitions and results. For simplicity, we consider 
throughout this section that q — (to, Wtg) G K x SU(n) is 
fixed and that Wq'. K — > SU(n) denotes the solution of the 
auxiliar system ([S])-© with initial condition Wq{to) — Wtg- 

Definition 2: [11] A point W £ Af" is called a limit 
point of Wq if there exists a real sequence {tm} such that 
limm^co tm = 00 and limm^oo Wq{tm) = W. The set of 
all limit points of Wq is called the limit set of Wq and is 
denoted by fl{Wq). 

Remark 3: Since SU(n) is a compact subset of Af", it is 
clear that il{Wq) is a non-empty subset of SU(n). 

Proposition 1: [11] limt^oo d{Wq{t), n(Wq)) = 0. 

The next 2 lemmas are essential in the proof of the 
important convergence result of Theorem |3] given below, 
which was inspired in the periodic version of LaSalle's 
invariance principle presented in [11] and in the ad-condition 
stabilization method of [6]. 

Lemma 1: Let W: M Af" be a continuously differen- 
tiable mapping such that limt^oo W{t) = 0. Suppose that 
{t„i} is a real sequence such that limm^oo = 00 and 
lim„j^oo W{tm) ~ W . Then, for every e G M, we have that 
lim^^oo W{t„^ + e) = W. 

Proof: Let e > and m e N. We have that 
W{tm + e) - W{Un) = /t*"^' W{t) dt. Thus, the inequality 
\\W{t^ + e) - W{tm)\\ < |e|sup,g[,^^,^+,] \\W{t)\\ holds. 
The assumptions then imply that lim,„^oo W{tm + e) = W . 
For e < 0, we can proceed in an analogous manner. ■ 



Lemma 2: Consider that W G ^{Wq), j G N and let 
1 < k < m. Assume that \m\t^ooWq{t) — and that 
\hnt^^\{Wq{t)Xl{t)Bl{t)X^) ^ 0, where is as in 
©. Then, \{WX^Bi{Q)X^) = 0. 

Proof: Let W G Vl{Wq). By definition, there exists 
a real sequence {t,,,.} such that lim,„^oo = 00 and 
\\m„i^oQWq{t„i) = W. Now, for each to G N, there 
exists £m G Z such that s„i — tm — (-mT G [0, T), where 
T > is the period of Xr and of i?^ (see Remark [Hi. 
Since [0,r] is compact, there exists a subsequence {sm;} 
in which limi^oo Sm; — € [0,r]. Let {tm,} be the 
corresponding subsequence of {tm}- Define the sequences 
{t^J and {s*^J as t*^^ ^ tm, - and s*^^ = Sm, - 0, 
respectively. We have that limt^oo Wq (t) = as well 
as limt^ooyiWqit)Xl{t)Bl{t)Xo^) = (assumptions). 
Therefore, by definition, limi_»oo Sm — 0, and Lemma [T] 
gives that lim,_,oo = limi^oo Wq(tm, -9) =W. 

Hence, the continuity and periodicity of Xr and of B-j, 
imply that lim.^o, V(M/,(t:;jXt(t;^Ji3^(t;^jXo,) = 
lim_^oo Y{Wq {t*m^ )Xl [s*^^ )Blis*^jXoo) 

YiwxlBiio)x^) ^0. m 

Theorem 3: Consider the subset E = {W G SU(n) : 
Y{WXI^bI{0)X^) = 0, for all j G N, 1 < fc < to}, 
where B^ is as in (O. Then, Mmt^ao d{Wq{t), E) = and 
E is non-empty. 

Proof: Due to Proposition [T] it suffices to prove that the 
non-empty limit set r2(Wg) of the solution Wq is contained in 
the set E. We remark that since V: A/" — ^ R is a continuous 
linear function, there exists c > such that |V(X)| < c||X||, 
for all X G Af". Furthermore, it follows from ©, (Ell-®, 
Remark [T] and the compactness of SU(n) that each of the 
mappings Xr, X}., Wq, Bl, Xr, X^, Wq, Bl is bounded, 
for every j G N, 1 < fc < to. 

Consider the functions a: M -> M, 6].: R x A/" ^ R, 
Pi'. R ^ R defined respectively as 

a{t) =Y{Wq{t)), foralKGR, 
Ix'^it, W) = YiWX^rit)Blit)X^), {t, W)eRx Af", 
/5fc (0 = K . Wq (t) ) , for all t G R, 

for j G N, 1 < fc < TO. We will prove by induction that 
lim (3i{t) = lim YiWq{t)Xlit)Biit)X^) = 0, (16) 

t — >oo t — >co 

for j G N, 1 < fc < TO. From (H))-® and the definition of 
feO, we have that Y{t,W) = ELi /l^fc(^' ^ ™d 
Yit,W) = 2j:Z^fMt,W)bl{t,W), where b°{t,W) = 
Y{WXUt)j:Zj!b1it,W)mBlit)X^)+blit,W) and 
fi,...Jm G K are non-zero, for {t,W) G R x Af". 
Since V is a non-negative function, we conclude that a is 
a non-decreasing function bounded from above such that a 
is bounded. Hence, liiUf^oo ct{t) exists and is finite. This 
relation along with Barbalat's Lemma (see e.g. [10]) give 
that limt^ooa(t) - Er=i ^9 W)' = 0. Thus, 
limt^oo/3g(i) = limt^ooY{Wq{t)Xl{t)B'j^{t)X^) = 0, 
for each 1 < fc < to, from which (O-® imply that 

lim Wq{t) ^ 0. (17) 



Now, consider the induction hypothesis 

lim (3l{t) - lim Y{Wg{t)Xl{t)Bi{t)X^) = 0, (18) 

t — >-c<; t — *oo 

for some j G N and a\\ 1 < k < m. We have that 
biit,W) = Y{WXl{t)j:T=if!b1it,W)HeBiit)X^) + 
bl'^^it, W), for all 1 < fc < m, {t, W) gRx Af". Straight- 
forward computations show that (31 is bounded because 
f3l{t) = bl{t,Wq{t)), for all 1 < fc < TO, t e M. Hence, 
( fT^ and Barbalat's Lemma imply that limt^ao l^l^^ (t) = 
\imt^^YiWq{t)X^,it)Bi+^it)X^) = 0, for 1 < fc < TO. 
We have thus proved that (fTSI l is true. At this moment, it 
is simple to prove that n{Wq) C E. Indeed, assume that 
W e n{Wg) C SV{n). Then, dlSll, dnli and Lemma|2]imply 
that V{WxI,bI{0)X^) = 0, for each j e N, 1 < fc < to. 

■ 

Lemma 3: Consider the subset F — {W G SU(n) : 
V{W) = for some e C such that |A^| = 

l,nr=i^' = l,^i(Ai) - ••• - 3(A„)}. Then, / e F and 
limt^oo d{Wgit),F) ^0. 

Proof: According to Theorem [3] it suffices to show 
the inclusion E C F. Let W E E C SU(n). It is a 
well-known result in linear algebra that W S SU(n) can 
be decomposed as W = Mdiag(Ai. . . . , A„)M'I', where M 
is unitary, Ai;...;A„ G C, \Xi\ — 1 and H^Li '^^ ~ 
Thus, V{W) = Eti^i^^) and Y{WXlBi{0)X^) = 
V(diag(Ai,...,A„)(XooM)tB^(0)(XooM)) = 0, for each 
1 < fc < TO, j G N. Since X^oM is unitary, it 
is clear that N: su{n) su(n) defined by N{Y) = 
{XaoM)^Y{XaoM), for every Y £ su(n), is a linear surjec- 
tive isomorphism. Now, by assumption, system ([T]i is regular 
and (O is satisfied. Hence, V(diag(Ai, . . . , A„)X) = 0, for 
every X S 5u{n), and thus V(diag(Ai, . . . , \n)D() = 0, for 
each 1 < ( < n, where Di = diag(z, ~i, 0, . . . , 0), D2 = 
diag(0,i, 0, ...,0),...,Aj-i = diag(0, -z) and 
Dn = diag(?, 0, . . . , 0, —i) are the canonical diagonal matri- 
ces of su(n). From the diagonal structure of Di, . . . , £)„, we 
conclude that Ai, . . . , A„ must satisfy 3(Ai) = • ■ • = 3(A„). 
This implies that W £ F and therefore E d F. ■ 

The proof of Theorem [T] is given below. 

Proof: It is clear that n = max (G) because I £ F. 
We will first show that G is finite. Let x £ G. Then, there 
exist Ai, . . . , A„ £ C such that x = X]j=i ^(-^j) with (i) 
n;=iA„ (ii) |A,| = 1 and (iii) 3(Ai) = ••• = 3(A„). 
Property (ii) implies that Aj = e'^^, for some 6j £ M, and it 
follows from (iii) that Aj = Ai = e'^i or A^ = e*('^-^i), for 
each 1 < j < n. Let rii be the number of j £ {1, . . . , n} 
such that \j = Ai and define ^2 = n — ni. Therefore, 
X — ni cos(0i) + ri2 cos(7r — 9i) = (rii — ^2) cos((?i) with 
f^i, (f^2 + 1) G {1) ...,«} and ni +712 = n. If rii = 712, then 
X = 0. Thus, assume that ni ^ n^. From property (i) we 
obtain that e»"i»ie*"^(''"''i) = 1. Hence, there exists fc e Z 
such that +^2(71 — = 2fc7r. This relation implies that 
B\ = (2fc — n2)'K I (n — 2n2). Note that ni,n2,k depend on 
X £ G and that n2, ni — 712 can only assume a finite number 
of values. If we show that 9i can only assume a finite number 
of values, we will have shown that the same holds for x £ G, 



which implies that G is finite. It is clear that the function 
77: Z ^ R defined as 7^{t) = cos((2^ - 712)71/(77. - 2712)), 
for all £ G Z, has period \n — 2n2\ > 0. Thus, the values 
assumed by 9i must be finite in number. 

Now, the convergence result will be shown. Recall that, 
for all X £ SU(7i), we have that -77 < V{X) < n and that 
V{X) = n if and only if X ^ I. Let 5 = max(G' \ {77}). 
Since G is finite, we have that 

(5 < a; < 77 => X = 77, for all a; e G. (19) 

Suppose that V{Wta) > S. Since SU(7i) is compact and 
V: A/" ^ M is continuous, it follows that V^|SU(7i) is 
uniformly continuous. Define the function a: M M by 
a{t) = Y(Wq{t)), for t £ R. Recall that, by construction, 
we have that Y{t,W) = T,T=iak{t,W)^ > 0, for all 
{t, W) £ R X M". Note that Uk and V are smooth, for 
each 1 < fc < TO. Since V is a non-negative function, 
we conclude that a is a smooth non-decreasing function. 
Therefore, ViWg{t)) > V{Wt„) > 6, for all t > to. The 
uniform continuity of V^|SU(77) then implies that there exists 
/i > such that 

\\X~Wg{t)\\ < /i V{X) > (5, for t>to, X £ Sl]{n) 

(20) 

(indeed, choose e = V{Wto) — 6 > 0). The convergence 
result of Lemma |3] means that 

Ve > 3T e M Vt > T 3a{t) £ F s.t. \\a{t) - Wg{t)\\ < e. 

Let e > and define e = min(e, /i). Thus, 

\ft>T 3a{t) £ F s.t. \\a{t) - Wq{t)\\ <e< fi, 

for some T £ R. Define T = max(T,to) and let t > f. 
Since a{t) £ F £ S\J{n) and V{F) C G, (EOl gives that 
5 < V{a{t)) £ G. However, 6 < V{a{t)) < n. Therefore, 
from il9[ . we obtain that V{a{t)) = n, which implies that 
a{t) — I. We have thus shown that limf^oo Wq{t) — I. ■ 
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